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ABSTRACT 

The different electron distributions in the hard and soft spectral states of black-hole binaries 
could be caused by kinetic processes and changing because of varying physical conditions in 
the corona. In presence of a magnetic field in the corona, the electron distribution can appear 
thermal, even when acceleration mechanisms would produce non thermal distributions. This 
is due to fast and efficient thermalization through synchrotron self-absorption. In this paper, 
we have analyzed data from 6 years of observations of Cygnus X-1 with the INTEGRAL ob- 
i-^ ■ servatory and produced 12 high-quality, stacked broad-band hard X-ray spectra representative 

of the whole range of spectral shapes observed in this source. We then fit these spectra with 
hybrid thermal/non-thermal Comptonization models and study the evolution of the physical 
parameters of the accretion flow across the spectral transition. In particular, we use the BELM 
model to constrain the magnetic field in the corona through its effects on the coronal emission. 
Indeed, the hot electrons of the X-ray corona produce soft (optical-UV) synchrotron radiation 
which is then Comptonized and may affect the temperature of the electrons (and thus the slope 
of the X-ray spectrum) through Compton cooling. We find that in the softer states, the emis- 
sion is dominated by Comptonization of the disc photons and the magnetic field is at most 
of the order of 10^ G. In the harder states, the data are consistent with a pure synchrotron 
self-Compton model, although a significant contribution of Comptonization of disc photons 
may not be excluded. If the non-thermal excess observed above a few hundred keV in the 
hard state is produced in the same region as the bulk of the thermal Comptonization, we ob- 
tain an upper limit on the coronal magnetic field of about 10^ G. If, on the other hand, the 
non-thermal excess is produced in a different location (such as the jet or a different part of the 
corona), the constraints on the magnetic field in the hard state are somewhat relaxed and the 
upper limit rises to ^ 10^ G. We discuss these constraints in the context of current accretion 
flow models. 

Key words: Gamma-rays: observations - accretion, accretion discs - black hole physics - 
radiation mechanisms: non-thermal - X-rays: binaries - stars: individual: Cygnus X-I 



1 INTRODUCTION 

Black hole binaries (BHBs), and Cyg X-1 in particular, show 
dramatic spectral variability. The most dramatic variations ob- 
served in BHBs are spectral transitions which occur, for most 
of these sources, on the time scales from weeks to days 
dZdziarski et al. 20021 ). In the simplest classification, BHBs, are 
observed in two main spectral states: in the hard state (HS) 
the source emits essentially radiation peaking around 100 keV 
( IGierlinski et al. 1997irZdziarski & Gierkinski 20041 while in the 
soft state (SS) the X-ray spectrum is a prominent ~1 keV black- 
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body plus a power-law tail (Zdziarski 2000 and ref. therein). The 
HS spectrum is believed to originate from thermal Comptonization 
in a hot electrons dominated cloud (the so-called "corona"; Eardley, 
Lightman & Shapiro 1975; Sunyaev & Titarchuk 1980). 

The black-body dominated spectrum observed in SS 
was associated with the optically thick accretion disc 
( [Shakura & Sunyaev 1973[ >, while the origin of the weak steep 
power-law component is still debated Cdziarski 2000t . In addition 
spectra with properties that are intermediate between that of the SS 
and HS are observed when the sources switches from one of the 
two stable state to the other. Terms such as "hard intermediate", or 
"soft intermediate" are often used to describe different flavors of 
those intermediate state l IBelloni et al. 20051 1 but there is no unique 
classification (see McClintock & Remillard 2006). 

The spectra of BHBs consist in several components. In the 
HS of Cyg X-1, the electron temperature and the Thomson op- 
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tical depth of the Comptonizing plasma were found to be typ- 
ically fcTo = 50-100 keV and r = 1-3 dGierlinski et al. 19971 
IFrontera et al. 200T| l. Although the scattering electrons have a pre- 
dominantly Maxwellian energy distribution in the HS, there are 
indications that the electron distribution may have some high- 
energy tail i.e., it is hybrid, thermal/non-thermal. The first evi- 
dence has been provided by CGiJO/COMPTEL observations at 
E > 500 keV of Cyg X-1 in hard state dMcConnell et al. 19941 
IMcConnell et al. 2002) . A slight indication was also inferred by 
CGRO/OSSE observations of GX 339^ ( Wardzihs kTet al. 2002b . 
Then, with the INTEGRAL mission, this feature has been ob- 
served many times during hard and/or intermediate/hard states 
in a number of BHBs and BH candidates (BHC; see Del 
Santo 2012 for a review on this topic), such as IE 1740.7- 
2942 dBouchet et al. 20091 1, GX 339-4 dPel Santo et al. 20081 
IJoinet et al. 20071 IDroulans et al. 2010t , Cyg X-1 (Malzac et al. 
2006; Cadolle Bel et al. 2006; Jourdain, Roques & Malzac 2012; 
Zdziarski, Lubinski & Sikora 2012). 

Alternatively, it is possible that the non-thermal excess and 
thermal Comptonization component are not produced by the same 
electron distribution but in separate regions of the corona with dif- 
ferent physical properties (Malzac 2012). Also, radio observations 
of Cygnus X-1 have resolved a compact relativistic jet (Stirling et 
al. 2001) and it was suggested that the high energy excess may be 
a component arising from it (Rahoui et al. 2011). This idea has 
been also proposed following the hard X-ray polarization measure- 
ments by Laurent et al. (201 1) obtained with the IBIS telescope in 
Compton mode. These authors found that, while the 250-400 keV 
spectrum is consistent with emission dominated by Compton scat- 
tering by thermal electrons and are weakly polarized (Pf < 20%), 
the second spectral component seen in the 400 keV-2 MeV band 
is strongly polarized (Pf = 67± 30%). They have argued that the 
MeV excess is likely to be produced in jet through synchrotron 
emission in a very coherent magnetic field. Such a high degree of 
polarization is indeed difficult to achieve through inverse Compton 
emission in the corona. Recently Jourdain et al. (2012) reported on 
a similar result by using the SPI telescope. They found that above 
230 keV the Cyg X-1 emission is indeed polarized, with a mean 
polarization fraction of 76± 15%. 

The shape of the MeV excess however, appears difficult to re- 
produce in the jet scenario (Zdziarski, Lubinski & Sikora 2012). 
Such a scenario would also imply that the MeV excess in the HS 
and the non-thermal power-law emission observed in the SS have 
a completely different nature since the radio jet is quenched in the 
SS. We note that although it is very plausible that the non-thermal 
excess originates in the jet, it appears very unlikely that the jet emis- 
sion could dominate the entire high energy spectrum of Cyg X-1 
(see e.g. Malzac, Belmont & Fabian 2009 and references therein). 

The X/7-ray SS spectrum of Cyg X-1 has been studied ex- 
tensively by simultaneous observations with ASCA, RXTE, Bep- 
poSAX, and CGRO during summer 1996 dPi Salvo et al. 20011 
IFrontera et al. 20 01). In addition to the dominating black-body, a 
long power-law like tail extending up to 10 MeV was discov- 
ered dMcConnell et al. 2002) . The high energy spectrum is well de- 
scribed by single Compton scattering off electrons having a nearly 
power-law distribution (Gierlinski et al. 1999, hereafter G99; Fron- 
tera et al. 2001). 

The different spectral states are usually understood in terms 
of changes in the geometry of the accretion flow. According to a 
popular scenario (see e.g. Done, Gierlinski & Kubota 2007 for an 
exhaustive review), in the SS a geometrically thin accretion disc 
(the standard Shakura-Sunyaev disc) extends down to the last sta- 



ble orbit and is responsible for the dominant black-body emission. 
This disc is the source of soft seed photons for Comptonization in 
small active coronal regions located above and below the disc. The 
magnetic field lines rise above the accretion disc through magnetic 
buoyancy, transporting a significant fraction of the accretion power 
into the corona where it is then dissipated through magnetic recon- 
nection dGaleev et al. 1979] l. This leads to particle acceleration in 
the corona. A population of high-energy electrons is formed which 
then cools down by up-scattering the soft photons emitted by the 
disc. This produces the high-energy non-thermal emission (see e.g. 
G99; Zdziarski et al. 2002). In the HS, the standard geometrically 
thin disc is truncated at distances ranging from a few tens to a few 
thousand gravitational radii from the black hole. In its inner parts, 
the accretion flow takes the form of a hot geometrically thick, op- 
tically thin disc (Esin, McClintock & Narayan 1997). In this hot 
accretion flow, the electrons are predominantly heated by Coulomb 
interaction with a population of hot ions and cool down by Comp- 
tonizing their own synchrotron emission and/or soft photons from 
the accretion disc. 

The different electron distributions in HS and SS could be 
caused by kinetic processes and changing because of varying phys- 
ical conditions in the corona (Malzac & Belmont 2009, hereafter 
MB09; Poutanen & Vurm 2009, hereafter PV09). These authors 
showed that in presence of a magnetic field in the corona, the elec- 
tron distribution can appear thermal, even when acceleration mech- 
anisms would produce non thermal distributions. This is due to fast 
and efficient thermalization through synchrotron self-absorption as 
first pointed out by Ghisellini, Guilbert & Svensson (1988). MB09 
and PV09, presented a rough 'fit by eye' of the average CGRO data 
of Cyg X-1 with this synchrotron boiler model which already pro- 
vided estimates of the magnetic field in the hard state under the 
assumption that the non-thermal excess is produced by electrons in 
the same zone as the Maxwellian component. In this paper, we will 
use the same model to perform a statistical fit of INTEGRAL data 
for the whole range of observed spectral shapes and also consider- 
ing the possibility that the non-thermal excess in the hard state may 
have a different origin. 

As one of the brighter galactic hard X-ray source, Cyg X-1 is 
a prime target for the INTEGRAL mission jWinkler et al. 2003) . It 
was extensively observed during the Performance Verification (PV) 
Phase of the mission. The on-board instruments offer an unprece- 
dented simultaneous broad-band spectral coverage, ranging from 3 
keV to several MeV. Thus a large amount of observing time (open 
time and core programme for a total of about 7 Ms) have been ded- 
icated to this target by INTEGRAL. A number of studies of the in- 
termediate spectral states have been performed (Malzac et al. 2006; 
Cadolle Bel et al. 2006); an X-ray flare during the hard state ap- 
pears to be coincident with the TeV emission detected by MAGIC 
(Albert et al. 2007; Malzac et al. 2008). 

Recent papers have also provided a refined estimation of the 
distance (d=1.86 kpc; Reid et al. 201 1) and the black-hole mass of 
Cyg X-1 (14.8M0;Orosz et al. 2011). 

In this work, we present the long term behaviour and a spectral 
variability study of Cyg X- 1 using the whole INTEGRAL data base 
available until Spring 2009. Rather than following the evolution of 
the source chronologically, the aim was to extract high quality data 
representative of the source in a given 'spectral state'. For this pur- 
pose, we stacked all the pointings with similar hardness ratios to 
produce 12 averaged spectra spanning the whole range of observed 
spectral shapes. The details of how this was done are described 
in Section |2] These 12 broad band (3 keV-1 MeV) spectra were 
then fit with two different hybrid thermal Comptonization models 
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Figure 1. From top to bottom: ASM count rate in 1.5-3 keV, ASM hardness 
ratio (3-6 keV/1 .5-3 keV), IBIS/ISGRI count rate in 20-40 keV, IBIS/ISGRI 
hardness (40-80 keV/20-40 keV). 

namely EQPAIR ( |Coppi 1999l l and BELM (Belmont, Malzac & Mar- 
cowith 2008; hereafter BOS) that are described in Section |3] This 
allows us to constrain the physical conditions in the corona and 
determine how the physical parameters change during the spectral 
evolution. In particular, fits with BELM sets quantitative constraints 
on the strength of the magnetic field in the corona of Cygnus X- 
1 for the first time. These results are described and discussed in 
Section|4] 
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Figure 2. Number of the INTEGRAL pointings considered in this paper vs 
IBIS/ISGRI hardness ratio (40-80 keV/20-40 keV). 

For the purpose of spectral analysis, we have built our data-set 
by selecting all IBIS observations including Cyg X-1 in a partially 
coded FOV of 15° x 15° where the instrument response is well 
known. Based on this selection, the total number of IBIS and SPI 
science windows (SCWs) results as 1024 (Tab.[T) for a total expo- 
sure time of about 3.3 Ms (a single SCW can last from 800 up to 
4000 seconds). The total SPI exposure time is essentially the same 
as that of IBIS, even though a further analysis excluding SPI prob- 
lematic pointings (i.e. solar flares, radiation belts) leads to a slightly 
lower SPI effective exposure time than IBIS. Because of the smaller 
FOV of the X-ray telescope JEM-X, the related selected pointings 
(FOV of 3.5° X 3.5° has been considered) are 364 (see Tab.[B. In 
particular, 44 SCWs for JEM-X2, operating at the beginning of the 
mission and 320 for JEM-Xl operating since revolution 181. 

The IBIS scientific analysis, focused on the low energy de- 
tector ISGRI dLebrun et al. 2003 ). and the JEM-X analysis have 
been performed using the INTEGRAL off-line analysis soft- 
ware, OSA v.9 dGoldwurmetal. 20031 ICourvoisier et al. 20031 1. 
The IBIS/ISGRI light curves (Fig. [Hi have been extracted in two 
energy bands (20-40 keV, 40-80 keV). JEM-X and IBIS/ISGRI 
spectra, extracted SCW by SCW, extend from 3 to 20 keV and from 
20 to 500 keV, respectively. 

To account for change over time of the energy response, IBIS 
is calibrated in different periods using dedicated observations of 
the Crab. Then for each Crab observation an ancillary response file 
(ARE) is estimated and associated with a given interval of revolu- 
tions in the INTEGRAL data (see IBIS User Manual). Since spectra 
in our data analysis (see Section fT2\ often correspond to obser- 
vations collected in different periods, i.e. associated with different 
ARE files, we evaluated for each averaged spectrum an ARE file 
obtained by a weighted mean of the effective areas in the different 
periods involved. 



2 OBSERVATIONS AND DATA ANALYSIS 

We have analyzed six years of Cyg X-1 observations per- 
formed by INTEGRAL until April 2009 (revolution 794) with 
the three coded-mask telescopes IBIS dUbertini et al. 20(Ht , SPI 
dVedrenne et al. 20031 1 and JEM-X dLund et al. 2003) . We started 
from March 2003 (revolution 52), thus we excluded the perfor- 
mance verification phase observations (1 Ms) due to some peculiar 
instrument settings. Results on Cyg X-1 obtained during the cali- 
bration phase have been extensively reported (Bouchet et al. 2003; 
Pottschmidt et al. 2003; Bazzano et al. 2003). 



2.1 SPI data treatment 

The signal recorded by SPI camera is composed of contributions 
from sources in the FOV convolved with the instrument response 
function plus the background. A system of equations is to be solved 
to determine sources and background intensities. In order to reduce 
the number of uncertainties and hence to maximize the signal-to- 
noise ratio of each sources, a few assumptions are introduced. 

The background count rates of the 19 Ge detectors (unifor- 
mity maps) are assumed to be stable on timescale of ~6 months, 
while the global normalization factor is determined on ~6 hours 
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Table 1. Observations log of the twelve final data-set used for spectral mod- 
eling. The IBIS harness ratio intervals (HR), number of pointings (SCW) 
with IBIS and SPI and with JEM-X (within brackets), IBIS and SPI observ- 
ing exposures (Expl), JEM-X exposure (Exp2) and the power-law slope (F) 
of the 30-100 keV IBIS/ISGRI spectra are shown. 



HR 


SCW (JEM-X) 


Exp it 


Exp2t 


r 




number 


[ks] 


[ks] 




0.37-0.41 


7(1) 


15.4 


2.2 


2 6+0.04 
^■"-0.03 


0.43-0.48 


12(2) 


36.8 


6.5 


2.32 ±0.02 


0.48-0.50 


35 (9) 


111.5 


27.8 


2.29 ±0.01 


0.50-0.53 


59 (8) 


181.1 


22.9 


2.23 ±0.02 


0.53-0.56 


82 (11) 


236.8 


34.0 


2.17 ±0.02 


0.56-0.58 


38 (19) 


101.0 


51.1 


2.08 ± 0.02 


0.58-0.61 


113 (40) 


343.4 


113.7 


2.02 ± 0.02 


0.61-0.64 


97 (27) 


287.7 


80.3 


1.98 ±0.02 


0.64-0.66 


114(29) 


944.2 


34.7 


1.81 ±0.02 


0.66-0.68 


340 (72) 


781.4 


147.8 


1.79 ±0.02 


0.68-0.70 


157(16) 


402.0 


34.2 


1.74 ±0.02 


0.70-0.72 


9(4) 


26.1 


8.0 


1.73 ±0.02 



3.0 [ 



2.5 



f-H 2.0 



1.0 L 
0.3 



0.4 



0.5 0.6 
HR 



0.7 



0.8 



Figure 3. Congelation between the spectral photon index of the IBIS spectra 
fit with a simple power law in the range 30-100 keV and the hardness ratio 
40-80/20^0 keV. Dashed lines indicate the twelve data groups used for the 
final spectral analysis (see text). 



timescale. For the sources, the timescales are chosen in function 
of the source intensity and temporal behavior: the faint sources 
are considered as constant. More information and details on SPI 
data reduction can be found in Bouchet et al. (2008) and Bouchet, 
Roques & Jourdain (20 11). 

At energy around 1 MeV, high energy particles saturate the 
electronics and can generate false events. The trigger signal issued 
by the Pulse Shape Discriminators is used to select events between 
650 keV and 2.2 MeV as explained in Jourdain & Roques (2009). 

In order to maximize the signal-to-noise ratio, we used a cata- 
log containing sources detected above la in the whole time interval 
spawned by the data. The catalogue contains 10, 6, 5 and 1 sources 
in the 25-50 keV, 50-100 keV, 100-300 keV and 300-600 keV en- 
ergy bands, respectively. In addition Cyg X-1 is time- variable up 
to ~400 keV, while above ~400 keV, the signal-to-noise ratio is 
better when assuming Cyg X-1 as constant. Above 500 keV, we 
can derive a better uniformity map using high latitude exposures 
(telescope pointing axis satisfying \b\ > 30°) as they contain no 
significant source emission dBouchet et al. 2010t . Comparison of 
the fluxes obtained with this uniformity map and one that is deter- 
mined by the fitting algorithm (Bouchet et al. 201 1) shows that the 
difference in determined fluxes is well smaller than the error bars 
and hence fluxes are fully compatible. In this paper, we use SPI 
spectra from 22 keV up to 1 MeV (when possible because of the 
statistics). 

2.2 Long term behaviour and spectral variability 

As a first step in characterizing the spectral variability (i.e. spectral 
states) of Cyg X-1, we computed the hardness ratio (HR) defined 
as the count rate ratio 40-80 keV to 20-40 keV for each one of 
the 1024 IBIS pointings (Fig. [T] panel 4). Cyg X-1 was most of 
the time in HS so that observations in hard states are significantly 
more numerous than the ones in soft states (Fig. O. During the 
first three years of the INTEGRAL monitoring, Cyg X-1 showed 
a very variable activity in term of both flux variation and spectral 
transition. Since middle of 2006, the source entered in an almost 
steady hard state, characterized by flux variability not combined 
with simultaneous spectral variability (see the two flares around 
MJD~54000 and MJD~54600 and the corresponding flat HR in 



Fig. [T] panel 3 and 4, respectively). For comparison with the soft 
X-rays, we plot the RXTE/ ASM long term behaviour, as well as the 
hardness ratio (Fig.[T] panels 1 and 2). 

We have fitted individually all 1024 IBIS/ISGRI spectra from 
20 keV to 80 keV with a simple power law model and verified that 
at a SCW level IBIS spectra in this band are described by a power- 
law with a null hypothesis probability higher than 5%. Therefore, 
the hardness ratio we have defined can be used as a spectral state 
indicator. 

The 1024 IBIS/ISGRI by-pointing spectra have been first aver- 
aged within each HR bin (HR bin-size=0.01), resulting in 31 spec- 
tra covering the range HR=0.37-i-0.72. These spectra have been fit- 
ted in the range 30-100 keV with a simple power-law. The correla- 
tion photon index vs. HR is shown in Fig.[3l Then, for the second 
step in averaging our IBIS spectra, in order to obtain a better statis- 
tics at high energy (above 100 keV), we have used two combined 
criteria. We averaged data from neighboring HR bins with com- 
patible (90% confidence level) 30-100 keV power-law slopes (first 
criterium; see Fig. [3]l until the error on F is reduced down to at 
least 0.02 (second criterium; Tab. [T). At the end we have obtained 
12 IBIS spectra collected with distinct shape spanning most of the 
range of observed spectral diversity of Cyg X-1. 

The JEM-X spectra have been then averaged following the 
IBIS selection. We checked the JEM-X spectral variability inside 
our IBIS HR bins and found that spectra inside the HR bin we have 
chosen are compatible with a constant slope in the JEM-X band. 
The SPI spectra have been extracted from the same 12 IBIS SCW 
groups. 



3 MODELS 

The twelve averaged broad-band spectra have been fitted with 
XSPEC V. 11.3.2, using two Comptonization models, namely EQ- 
PAIR (Coppi 1999, G99) and BELM (BOS, MB09), plus disc reflec- 
tion and relativistic iron line emission. The best-fit parameters of 
the twelve averaged broad-band spectra are shown in Tab. [3](EQ- 
PAIR) and Tab.|4](BELM). All quoted errors are at 90% confidence 
level. 
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3.1 Coronal emission models 

EQPAIR and BELM model the emission from a homogeneous, 
isotropic, spherical cloud of fully ionized matter in interaction 
with the ambient radiation field. Not only do they compute the 
observed photon spectrum of the corona, but also they give the 
energy distributions of leptons in the corona. Those results are 
computed self-consistently according to the following microphysi- 
cal processes: lepton-lepton Coulomb collisions, Compton scatter- 
ing, photon-photon pair production, pair annihilation, self-absorbed 
bremsstrahlung radiation, and self-absorbed synchrotron radiation 
(the latter for BELM only). The leptons distributions are generally 
found to be composed of a low-energy thermal population and an 
additional high-energy non-thermal population, so that those kind 
of models are usually called hybrid models. The main physical pa- 
rameters that characterize such corona are its typical size, density, 
magnetic field, and how it is energized. However, the spectral prop- 
erties depend more explicitly on non-dimensional combinations of 
these quantities. 

For instance, the dependence on the absolute coronal size is 
negligible for most parameter regimes, and the radius is set to i? = 
5 X 10^ cm, which corresponds to a region of about R ~ 20-Rg for 
a 15M0 black hole. 

The proton density Up is represented by the proton Thomson 
optical depth Tp = UpGTR, where ctt is the Thomson cross sec- 
tion. The Thomson optical depth Tp provided by the two models 
is related to the ionization electrons only. The total optical depth 
(Ttot) is the sum of the optical depth of e^e~ pairs plus Tp. 

The source luminosity L is described by the non-dimensional 
compactness parameter ( [Lightman & Zdziarski 1987] l: 



(1) 



where me is the electron mass, and c the speed of light. 

Photon escape can only be dealt with exactly by Monte-Carlo 
or complete radiation transfer simulation. Instead, it is modeled 
here through an approximate escape probability formalism (e.g. 
Lightman, Zdziarski & Rees 1987; Lightman & Zdziarski 1987; 
Stem et al. 1995). EQPAIR uses an escape probability based on 
pure photon scattering, whereas the escape probability used in the 
present version of BELM is based on both photon scattering and ab- 
sorption (Sobolev 1974; VP09), which can lead to significant dif- 
ferences for low energy photons and large optical depths. Although 
results with the two models can give slightly different results for a 
given set of parameters, the spectral dependence on parameters is 
identical in both models. 

To compensate for losses associated to photon escape, energy 
can be injected in the corona (i) directly as external photons from 
the accretion disc, or indirectly to the leptons populations through 
two channels: (ii) thermal heating or (iii) non-thermal acceleration. 

(i) External photons are assumed to originate from a optically 
thick, multi-temperature accretion disk extending from the mini- 
mal stable orbit of a non-rotating black hole J?in = (SRg out to 
Rout ~ 10^ J?g. The emitted multi-color black-body spectrum is 
then essentially determined by the temperature Tin at the inner ra- 
dius. As for the source luminosity, the power injected as soft pho- 
tons is characterized by the soft compactness: 1^ (also named Zbb in 
EQPAIR) . 

(ii) Thermal heating provides energy to the pool of thermal lep- 
tons with a compactness parameter /th- 

(iii) Non thermal processes such as magnetic reconnection 
and/or shock acceleration are expected to generate power-law like 



distributions of particles. This is modeled as an injection of leptons 
with a power-law distribution /(7) oc characterized by the 

minimal and maximal Lorentz factor 7min = 1-3 and ^max = 10^ 
between which the acceleration takes place, the power-law index 
Finj, and the total injected power described by the compactness 

^nth- 

= ^th + ^nth corresponds to the total power supplied directly to 
the plasma. In steady state, all injected power goes out as photons 
and I = hh + ^ith- As far as the compactness I does not reach 
too large values {I < 100), the spectral shape is rather insensitive to 
the total compactness and to the absolute values of Is, Ith, and Inth- 
However, it depends on the channel used to injected this power, i.e. 
on ratios such as Ih/ls and Inth/lh- For this reason, it is customary 
to fix Is to some reasonable value when fitting data with EQPAIR. 

The main difference between these two models is that in ad- 
dition to all physical process mentioned earlier, the BELM model 
allows for self-absorbed synchrotron radiation resulting from the 
gyro-motion of leptons around a tangled magnetic field, B (BOS). 
The magnetic field intensity is set through the magnetic compact- 
ness parameter: 



^Rf 



(2) 



Since the shape of the spectrum depends on the ratio Ib/1 
rather than on the absolute value of Ib we computed the BELM 
code leaving Ib/Iii as free spectral parameter. 

Also, in the BELM code the thermalization process is treated 
self-consistently while in EQPAIR the electron distribution is as- 
sumed to be purely thermal at low energy (see MB09 for details). 

The computation time with BELM is comparatively longer 
than with EQPAIR and the former cannot be used for inline fitting. 
Instead, a table of spectra was computed and used in XSPEC as a 
table model. 



3.2 Reflection on the accretion disc 

Compton reflection has been modeled by using the viewing-angle- 
dependent Green's functions approximation obtained by Magdziarz 
& Zdziarski (1995) for an isotropic point source (or, equivalently an 
optically thin corona) above a slab. We treat Q/2n as a free param- 
eter, where is an effective (i.e., corresponding to the observed 
strength of reflection) solid angle subtended by the reflector. The 
reflection is accompanied by an Fe fluorescence Kq line centered 
at an energy, -Eiinc- 

When the reflection comes from a fast rotating disc in a strong 
gravitational potential, Doppler and gravitational shifts become im- 
portant. These effects are usually approximated by convolving both 
the reflected component and the Fe Ka line with the Schwarzschild 
disc line profile of Fabian et al. (1989; DISKLINE in XSPEC). 

In calculating the relativistic distortion, we consider a range 
of radii from i?in=6 to J?out=1000 gravitational radii . The disc re- 
flection component is calculated for neutral material with standard 
abundances. We assumed an inclination angle of 50°. 

Due to the energy resolution and line sensitivity, the iron line 
study cannot be performed with JEM-X. The line energy was there- 
fore imposed at 6.7 keV. The inclination angle, inner radius of the 
disc and the disc emissivity law were fixed at the same values as in 
the Compton reflection model. The normalization is therefore the 
only free parameter for the iron line emission. 
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Figure 4. EQPAIR spectral parameters reported in Tab. [3]vs V reported in 
Tab. [T] In the bottom panel fluxes of the two principal components, disc 
black-body (squares) and Comptonization (triangles), are plotted. 



Figure 5. Joint JEM-X, IBIS and SPI energy spectra of Cyg X-1 during four 
different spectral states fitted with the thermal/non-thermal hybrid Comp- 
tonization model EQPAIR plus DISKLINE (see Tab.|3). 



4 RESULTS 

4.1 Unmagnetized model 

As a first study, we investigate here a model with no mag- 
netic field (Zs=0), so that we performed the analysis with EQ- 
PAIR. According with previous studies (Zdziarsk i et al. 20021 
IWilms et al. 2006t , we also have found that the bolometric lumi- 
nosity of Cyg X-1 does not change significantly from SS to HS. 
Therefore, unless the size of the corona changes dramatically dur- 
ing the transition, the total compactness is not expected to vary 
much ( [Poutanen & Coppi 1998) . However the distribution of the 
radiation changes dramatically. In the hard state most of the lumi- 
nosity comes from the corona while in SS most of the radiation 
is in the black-body disc emission. We therefore fixed Zs = 10 
when fitting soft and intermediate/soft states (States \-4 in Tab.[3j 
G99) and Is = 1 when fitting the intermediate/hard and hard states 
(States 5-12; Ibragimov et al. 2005). 

The best fit parameters are shown in Tab.|3] We note that spec- 
trum 1 shows spectral parameters of the canonical soft state, ac- 
cording also with the hard X-ray power-law slope (Tab. [T), while 
spectra 9-12 show typical spectral parameters of the hard state (see 
Zdziarski 2000; Zdziarski & Gierlinski 2004). Spectra 2-8 may be 
defined as intermediate. In addition, we call intermediate/soft all 
spectra where a disc black-body is required by the fitting procedure 
and intermediate/hard the spectra where it is not possible to give 
constraints to the inner disc temperature. 

In Fig. |4] we show the relation between the best-fit model pa- 
rameters and the hard X-ray power-law slopes. Four significantly 
different energy-spectra of Cyg X-1 are plotted in Fig. [5] while in 
Fig. [7] we show the total best-fit models for each of the twelve 
selected spectra. 

In agreement with previous studies (G99) our fits show that 
the electron-positron pairs constitute only a small, if not negligible, 
fraction of the Comptonising leptons. An interesting trend (shown 
in Fig.|4j is that the Thomson optical depth is lower in softest states 
(Tp ~ 0.2 in SS). It increases significantly from SS (spectrum 1) 
to intermediate states (spectra 2-8) and is approximately constant 
in the hard states (rp ~ 1.4, spectra 9-12). This change by a factor 
of 7 could be associated either to a change in electron density or 
in the size of the corona (since r oc UeR). If this is the size of 
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the corona that changes then this would imply that the corona is 
larger in the HS than in the SS by a factor of ~7. This would be in 
agreement with the decrease of the inner radius of the accretion disc 
and faster rapid variability in the SS which both suggest a smaller 
emitting region in the SS. But the electron density could also drive 
the change in Thomson depth. This could be the case if the material 
in the corona gradually condensates into the accretion disc when 
the source evolves toward softer states. Another possibility would 
be that the coronal material is evacuated during the ejection event 
(optically thin radio flare) at the spectral transition (as suggested 
by Rodriguez & Prat 2010). However we see on Fig. |4] that the 
evolution of Tp is gradual. It starts in the intermediate/hard states 
and continues in the soft state, this appears in contradiction with 
a scenario in which most of the corona would be suddenly ejected 
at the transition. If the corona is ejected this has to be a gradual 
process, that might be associated with the presence of strong disc 
winds evidenced in the SS of X-ray binaries (see e.g. Ponti et al. 
2012). 

We find that the inner disc temperature fcTin decreases from 
0.6 and 0.4 keV between spectrum 1 and 4 (soft state to interme- 
diate/soft state, see Tab.O. This is combined with a decrease by a 
factor of four of the flux of the measured thermal disc component. 
For harder spectra, the disc temperature is so low that it cannot be 
constrained with our data. We therefore imposed an inner disc tem- 
perature of 0.3 keV when in intermediate/hard states (spectra 5-8) 
and 0.1 keV when in canonical HS (spectra 9-12). The temperature 
of the corona has the opposite behaviour and is larger for harder 
spectra, except for spectrum 1. This change of coronal temperature 
is driven by a change in the ratio of the electron heating rate to 
soft cooling photon flux in the corona (oc Ih/h) which is larger at 
small F. The fact that the temperature appears higher in the softer 
spectrum 1, is related to the significantly smaller optical depth in 
soft state. In fact, even if the temperature is higher in spectrum I, 
the Compton parameter y oc ttTe which is expected to vary ap- 
proximately like Ih/ls is the lowest. Then computing the absolute 
disc and Comptonized fluxes (respectively Ftb and Fcompt, see 
Tab. [3} can provide indications on whether this evolution is driven 
by changes in the flux of soft cooling photons from the disc (oc Is) 
or changes in the electron heating rate (oc Ih). We find that the ther- 
mal disc flux changes by more than one order of magnitude while 
the Comptonized radiation decreases by almost a factor of two. The 
softening we observe in the states 1-4 appears to be caused by a dra- 
matic increase in the disc thermal flux in the corona associated with 
a modest reduction of the electron heating rate. 

Non-thermal acceleration appears to be required to fit most 
spectra. When trying to fit with a non-thermal compactness frac- 
tion fixed to /nth//h = 0, only for the spectrum 12 we obtained an 
acceptable xt (i-e. 0.98). However, for this last spectrum the statis- 
tics is very low and residuals show that an additional non-thermal 
component is required at high energy. For the other spectra, sig- 
nificant Inth/h-i fractions are required. This indicates that a non- 
thermal electrons population may also be present in the corona (or 
elsewhere), together with the thermal population even during the 
hard states. This is in agreement with previous studies performed 
on intermediate and hard states of BH binaries, as also reported in 
Sec.[T] In soft and intermediate/sost states (1-4) though, it was not 
possible to constrain Inth/lh and this parameter was fixed at 0.99. 

From our spectral fits, from the softest (spectrum 1) to the 
harder states (up to spectrum 9) Finj parameter changes from 3.5 
to 2.0. We are not able to give constraints on this parameter for 
spectra 10-12, so that we fixed it to 2. Such a change in the slope of 
the injected non-thermal electrons suggests a different acceleration 
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Figure 7. Total models obtained by fitting all twelve broad-band spectra 
with EQPAIR plus DISKLINE (see Tab.|3) 



Table 2. Fitting parameters of the BELM table model used in XSPEC and 
their properties (minimal and maximal values, number of interpolation val- 
ues, and grid style). Other parameters are set to: I = 10, Is = I — /nth. 
^th = 0, 7mi„ = 1.3, 'ymax = 10^, i? = 5 X 10^ cm. 
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mechanism in SS and HS and may support the view that the non- 
thermal excess in the HS has a different origin (e.g. the jet). We note 
however that the results of MC02 suggest a different evolution (i.e. 
Finj larger in the HS). It is possible that the soft injection we infer 
in the SS is related to the lack of photon statistic above 100 keV 
and to the fact that the non-thermal compactnes fraction 'nth/'h 
was fixed when fitting the softest spectra. 

We checked the need for the reflection component freezing 
Q.I2-K at 0. The results show that the are consistently worst (F- 
test probabilities < lO"**), confirming that introducing the reflec- 
tion component improves significantly the fits. Our fits suggest that 
the reflection component is significantly higher only in the softest 
state (spectrum 1) of our data set, while it is almost constant within 
the errors in the other states. This is a known result which can be in- 
terpreted as follows. In the HS the reflecting disc intercepts a small 
fraction of the coronal radiation since the disc is truncated at a large 
distance from the black hole; in the SS the accretion disc extends 
down to the last stable orbit and is sandwiched by the illuminating 
corona (Done et al. 2007). 

4.2 Magnetized, one-zone model with pure non-thermal 
acceleration 

If the corona is magnetized, it must emit synchrotron emission that 
can add up to disc photons and be Comptonized. In this section 
we aim at answering the question: how large can be the magnetic 
field before synchrotron emission makes the modeled spectra in- 
compatible with the data? It was shown in MB09 and PV09 that 
processes such as the synchrotron boiler and Coulomb collisions 



can thermalize particle distribution efficiently, even when particles 
are accelerated trough pure non-thermal mechanisms. To minimize 
the number of parameters, we focus here on such a model with pure 
non-thermal acceleration. The alternative case of thermal heating 
will be discussed in the next subsection. 

Using BELM, we have computed a table of more than 5 x 10^ 
spectra and 5 fitting parameters: Inth/l, Ib /'nth, Tp, Tinj, and kTm 
(see table|2j- This pure SSC model does not include any thermal or 
Coulomb heating (hh = Ic = 0). As shown in MB09, the spectrum 
shape depends on ratios such as Zs/'nth or Inth/l and not on the 
total compactness of the source. Therefore the best fit values for 
these ratios are also insensitive to the exact value of I. The total 
compactness was thus set to the typical value I = 10 and Is was 
computed such as Is + Znth — I = 10. 

On the contrary, the magnetic field intensity depends on the 
exact amplitude of the magnetic compactness Ib- Once the spectra 
were fitted and the flux derived from XSPEC, the real source com- 
pactness 'obs was computed using Eq. 1 and assuming a distance 
of d = 1.86 kpc. The exact magnetic compactness was then scaled 
linearly to l^hs as: 



where I b /'nth and Intb/l are provided by the fit and insensitive to I. 
The magnetic field B was then estimated from equation assum- 
ing a size of the corona as R = 20Rg . Alternatively, the observed 
magnetic field can be compared to the value of magnetic field in 
equipartition with the radiative energy density. In the limit of low- 
energy photons (hu < mc(?), the corresponding compactness ratio 
is (Lightman & Zdziarski 1987): 

Ib _ Ib 'nth 47r/3 
Zsr Inth I I + tt/S 

This ratio is again insensitive to the choice for I and is also inde- 
pendent of uncertainties on the source size and distance. 

Tab. |4] shows the best-fit parameters. Fig. [8] presents a selec- 
tion of four different INTEGRAL spectra fitted with BELM, while 
Fi gs . 1 1 1 and [9l displav respectively the photon spectra and the par- 
ticle distribution, for the corresponding models. Overall this model 
provide good fit to the data (Tab.|4l(. This confirms that the data are 
compatible with pure non-thermal acceleration models. 

The parameters Tp, fcTin, Znth (or Ih/ls) and the reflection am- 
plitude have a similar behaviour as in the fits with the unmagnetized 
model EQAIR, for similar reasons (Fig.[6]l. The fact that Inth/l is 
lower in softer state is equivalent to the trend observed for the ratio 
Ih/ls- We note that when using BELM the electron acceleration in- 
dex Finj is better constrained by the data in the magnetized model 
since it also affects the number of synchrotron emitting electrons 
which in turns control the temperature of the Maxwellian electrons 
and the slope of the X-ray spectrum. The trend for a softer injec- 
tion index in softer state is less clear in the fits with this magnetized 
model. 

From the results of our fits, it is clear that in the soft state, 
the accretion disc photons represent the main source of soft seed 
photons (see Fig. llOt. In the harder states (spectra 6-12), it is not 
possible to obtain any constraints and the synchrotron photons may 
dominate. We have set Is — when fitting these spectra and ob- 
tained good fits. Based on the observation of correlations between 
X-ray flux and spectral index in several black hole binaries in the 
HS, Sobolewska et al. (2011) have argued that synchrotron seeds 
should dominates at luminosities below about a percent of the Ed- 
dington limit, while disc seed photons should dominate at higher 
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Figure 8. Joint JEM-X, IBIS and SPI energy spectra of Cyg X-1 during four 
different spectral states fitted with the BELM model with pure non-thermal 
acceleration plus DISKLINE and REFLECT (see Tab.|4). 

luminosities. It is interesting to note that Cyg X-1 has a stable lu- 
minosity located around this transition. It is therefore likely that 
both the accretion disc and the magnetic field contribute to seed the 
Comptonization process. 

Nevertheless, since it is not possible to constrain their respec- 
tive contribution from our data, assuming the extreme case of pure 
synchrotron seed photons allows us to set an upper limit on the 
magnetic compactness. Indeed the total flux of seed photons (syn- 
chrotron-l-disc) determines the slope of the X-ray spectrum (via en- 
ergy balance of the thermal electrons). Since the flux of synchrotron 
photons increases with I b , and the spectral slope is well constrained 
by the data, assuming that there are no disc photons is going to max- 
imize Ib- In the softer states the spectra are consistent with Ib = 
(since they can be fit with an unmagnetized model) but again we 
are able to obtain an upper limit on Ib- 

We note that the jump in Ib /Inth and B between spectra 5 
and 6 is due to the fact that we had to turn off the disc seeds pho- 
ton when fitting spectra 6-12. This makes the upper limit on Ib 
less restrictive. The presence of this jump suggest that disc photons 
continue to play a role even in the harder spectra. 

From Tab.El we see that for a corona of size R = 2QRg , these 
upper limits range from ~ 10^ G in the hardest HS to ~ 3 x 10^ G 
in the soft state (see also Fig.jSJ. This implies strongly subequipar- 
tition magnetic fields (Ib/Ibh < 1) in the spectra intermediate 
and hard. As discussed in MB09 and Droulans et al. (2010) , this 
has important consequences for dynamical models of the corona. 
Namely, this rules out any model based on a X-ray corona powered 
through magnetic dissipation in the hard state and this may also be 
a problem to drive the powerful jet that is observed. 

4.3 Magnetized, two-zone model with pure thermal heating 

In the previous one-zone model, we obtained upper limits on the 
magnetic field intensity assuming Synchrotron Self-Compton emis- 
sion (SSC) and pure non-thermal acceleration. This was moti- 
vated by the overall non-thermal nature of the softer spectra and 
by the high energy emission above the cutoff in the harder states 
(HR>0.64), none of which is reproduced with purely thermal lep- 
ton populations. In the hard state, these constraints were mostly 
governed by the non- thermal, highest energy data points. How- 




(7-1) m^c^ (keV) 

Figure 9. Electron energy distributions obtained in the best fit models of 
Figs. [To] and [s] The different colours corresponds to those of the spectra 
shown in Fig. [8] 



ever, it has been suggested that the high energy emission in the 
hard states could originate from a different region of the accret- 
ing system, i.e. the jet dLaurent et al. 20lTl IZdziarski et al. 20T2t 
or another component of the corona (Malzac 2012). Then, the hard 
spectra below the cutoff could also be reproduced by SSC models 
with pure thermal heating. In this case the exact heating mechanism 
is of little importance for spectral modeling, because whatever the 
heating mechanism the resulting electron distribution is very close 
to a pure Maxwellian. In BELM this electron thermal heating mech- 
anism is assumed to be Coulomb collisions with a distribution of 
hot thermal protons (MB09). Then, the constraints on the magnetic 
field are expected to be different. 

To investigate this issue, we computed a new table model with 
pure Coulomb heating (Is = 0, /nth = 0, 1 = Ic = hh)- As before, 
the spectral shape does not depend on the absolute compactness of 
the system, and we set / = Zc = 1 to prevent degeneracy issues. 
Then we determined the proton temperature (kTp in Tab.O corre- 
sponding to this compactness and to the exact steady state electron 
distribution. The two free parameters of the table model are: Ib (20 
values between 0.1 and 10^), and r (20 values between 0.5 and 5). 
As for previous models, we also added a reflection component. We 
fitted the four harder spectra (HR>0.64; see Tab.[5j. 

As expected, the broadband spectrum can not be fitted with 
this pure thermal model, since it does not reproduce the highest 
energy data, which either requires a significant level of non-thermal 
acceleration or must originate in a different region. When excluding 
data above 200 keV, we obtained good fits with xi < 0.95. The 
results are shown in Tab.|5] The best fit Thomson optical depth and 
the reflection amplitude are comparable to those of the non-thermal 
model. 

Noticeably, the ratio Ib/Ic is larger by about 4 orders of 
magnitude in the thermal model than in the previous non-thermal 
model. Indeed, the spectrum below the cutoff is essentially gov- 
erned by the temperature of the thermal electron distribution. The 
cooling of thermal electrons is dominated by inverse Compton scat- 
tering on soft synchrotron photons. Non-thermal models produce 
stronger high energy tails in the electron distribution and thus tend 
to provide more synchrotron seed photons than thermal models for 
comparable I b ■ Hence the magnetic field has to be weaker in non- 
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Figure 10. BELM total models (black) corresponding to fits of the spectra shown in Fig.[8]and the components: Compton (red), pair annihilation/production 
(orange), synchi'otron (blue), soft black body photons (purple). Solid lines con'espond to positive contribution to the spectrum, dashed lines coirespond to 
negative contribution. 



thermal models in order to provide the same electron temperature 
and reproduce the data. 

Thus, by assuming that the highest energy emission originates 
from a different region, the constraints on the ratio Ib jlc estimated 
in this thermal model correspond to the most conservative upper 
limits. Following the procedure described in the previous subsec- 
tion, we obtain that the magnetic field in the hard states is stronger 
by a factor of 100 compared to that estimated in the SSC model 
(Tab. |5]l. The magnetic field is now allowed to be strongly super 
equipartition in the hard state with typical upper limits of the order 
of 10^ G. 

In the case of a pure Maxwellian plasma, the dependence of 
the Coulomb compactness (Zc) on electron and proton temperatures 
can be re- written as: 



(5) 



where (/(^c) = V79c and 9c = -^-^ is the electron temperature in 
units of the electron rest mass (see MB09 for details). Zc(obs) ~ 'obs 
is obtained from equation (TJ and the bolometric fluxes by the fits 
(Tab. |5}. We see that temperature of the ions is at most of a few 
MeV which is incompatible with standard two temperature accre- 
tion flow models (such as ADAF models). As discussed in MB09, 



this relatively low proton temperature is related to the large (i.e. 
> 1) Thomson optical depth required to fit the data. As discussed 
in Malzac (2012), standard hot disc solution do not produce such 
large Thomson depth (and low ion temperature). The magnetically 
dominated accretion flow model of Oda et al. (2010, 2011) may 
however produce such parameters. The strong magnetic field as- 
sumed in this model would be consistent with our upper limits only 
in the case of a pure thermal electron plasma with the MeV tail 
produced in the jet or in an other region of the corona. 



5 CONCLUSION 

We have used data from 6 years of observation of Cygnus X- 1 with 
INTEGRAL, to produce 12 high quality, stacked broad-band spec- 
tra spanning the whole range of spectral shapes observed in this 
source from the SS to HS. We have then fitted these spectra with 
hybrid thermal/non-thermal Comptonization models. We have used 
both a magnetized BELM and a unmagnetized EQPAIR models. For 
both models we found that the spectral changes at high energy are 
driven essentially by the strength of the soft seed cooling photons 
in the corona as the disc temperature (and inner radius of emission) 
changes. 
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During the spectral transition from HS to SS the power of the 
corona decreases by less than a factor of 2 while the disc emission 
increases at least by a factor of ~ 50. This appears to be associated 
with a gradual change of the Thomson depth of the corona which 
decreases from ~ 1.5 in HS to ~ 0.2 in SS, that may be caused 
either by a condensation of the corona into the disc or evaporation 
into an outflow, or compactification of the corona in the SS. 

For the first time we have performed a statistical fit of the 
Cyg X-1 data with the BELM model. This allowed us to test the 
effects of the magnetic field on the coronal emission. As explained 
in the Sec. 14.21 the estimate of B depends only on our choice for R 
and not on the assumed fixed value of / used in the fitting proce- 
dure. Within the one-zone model the uncertainty on R is by far the 
main source of error on our estimate of B, and the other parameters 
have none or very little effects. 

We found that in the softer states, the emission is dominated 
by Comptonization of the disc photons and the upper limit on the 
magnetic field is at most of the order of 3 x 10'' y^20Rg/R G. Such 
a magnetic field is energetically consistent with a magnetic dissi- 
pation process, such as reconnection, powering the corona. 

In the hard states the data are consistent with a pure Comp- 
tonized synchrotron emission model, although Comptonization of 
disc photons cannot be excluded. We consider that a mixture of 
both is the most likely possibility. However, fitting in the limiting 
case of SSC in the HS provides upper limits on the magnetic field 
intensity. 

If the non-thermal excess observed above a few hundred keV 
is produced in the same region as the bulk of the thermal Comp- 
tonization, the magnetic field is at most of a few 10* y^20Rg/R G. 
This is strongly subequipartion with radiation which is challenging 
for all dynamical accretion flow models. In this case, we also con- 
firm the qualitative results of MB09 and PV09 that the emission 
in all spectral states can be modeled with pure non-thermal elec- 
tron acceleration. Indeed, provided that the magnetic field is not 
too strong in the HS, the electrons thermalise at the observed tem- 
perature under the effects of Coulomb collisions and synchrotron 
self-absorption. Moreover, Wardzinski et al. (2002) found analyti- 
cally that if the high-energy tail observed in HS of BHBs is from 
the hybrid distribution, the magnetic field has to be strongly sube- 
quipartition, similar to our finding. 

If, on the other hand the non-thermal excess is produced in a 
different location, the constraints on the magnetic field in the hard 
state are somewhat relaxed and we obtain an upper limit of about 
2 X 10^^20i?9/i?G in the hard state. Such a large magnetic field 
would be consistent with the magnetically dominated hot accretion 
flow model of Oda et al. (2010) which also appear to produce the 
rather large Thomson depth (rp ~ 1.5) and low ion temperatures 
Tp ~ 10^" K, that we infer from the data in the hard state. 
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Table 3. The fit parameters of the joint JEM-X, IBIS and SPI spectra in twelve different spectral states. Fits have been peiformed with EQPAIR combined with DISKLINE. Values in parenthesis denote parameters 
fixed during the fits. One of the interesting fixed parameters is Zg (Zi-,b in EQPAIR) which is fixed at 10 for states 1-4 and at 1 for states 5-12. The inner disc temperature is fixed at 0.3 keV in the intermediate states, 
while it is 0.1 keV in the harder ones. See text for the parameters description. 
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Table 4. Same spectra than in Tab.[3]fitted with the BELM model with pure non-thermal acceleration plus DISKLINE and REFLECT. 
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Table 5. Spectra below 200 keV in HS fitted with a pure thermal heating model. kTp and feTp^oi^g) are the temperatures of the ions for Ic = 1 and Ic observed 
(Ic = /), respectively. 

HR Ib/Ic Tp Q/2TT x5(dof) FboI x IQ-* B Ib/Ibt kT^ kTp fcTp(„bs) 

[ergcm-^s-i] [Gauss] [keV] [keV] [MeV] 

0.64-0.66 68+3^ 1-83-0.02 0.12+J5 0.94(107) 3.5 1.9 x lO'^ 176.3 79.0 245 1.6 

0.66-0.68 AA+l 2.0 ±0.1 0.27to°4 0.94(58) 3.9 1.7 x lO'^ 110.3 70.2 219 1.5 

0.68-0.70 48^4^ 2.oloos 0.19lo°6 0.90(112) 3.0 1.5 x lO'^ 119.8 78.2 215 1.1 

0.70-0.72 64I2I 2.ll°-^ 0.14l°Jg 0.95(113) 4.1 2.0x10^ 157.1 72.3 183 1.3 
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